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Abstract 

We study the asymptotic behaviour of needlets-based approximate 
maximum likelihood estimators for the spectral parameters of Gaussian 
and isotropic spherical random fields. We prove consistency and asymp- 
totic Gaussianity, in the high-frequency limit, thus generalizing earlier 
results by Durastanti et al. (2011) based upon standard Fourier analysis 
on the sphere. The asymptotic results are then illustrated by an extensive 
Monte Carlo study. 
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1 Introduction 

In a recent paper, (see [T3], [2]) we investigated the asymptotic behaviour of a 
Whittle-like approximate maximum likelihood procedure for the estimation of 
the spectral parameters (e.g., the spectral index) of isotropic Gaussian random 
fields defined on the unit sphere S 2 . Under Gaussianity, it was indeed possible 
to establish consistency and a central limit theorem allowing for feasible infer- 
ence, under broad conditions on the behaviour of the angular power spectrum. 
These results, as many others concerning statistical inference on spherical ran- 
dom fields, have recently found strong motivations arising from applications, 
especially in a Cosmological framework (see for instance [39] and the references 
therein) . From the technical point of view, the asymptotic framework is rather 
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different from usual, as it is based on observations collected at higher and higher 
frequencies on a fixed-domain (the unit sphere) : even consistency of the angular 
power-spectrum estimator under such circumstances becomes non-trivial, and 
largely open for research, see also [55] . 

The procedure considered in [13] is based upon spherical harmonics and 
classical Fourier analysis on the sphere. Despite its appealing properties, it 
must be stressed that in many practical circumstances spherical harmonics may 
suffer serious drawbacks, due to their lack of localization in real space. This can 
make their implementation unfeasible, due to the presence of unobserved regions 
on the sphere (as it is commonly the case for Cosmological applications), and 
it may exclude the possibility of separate estimation on different hemispheres, 
as considered for instance by [BJ, [55]. In view of these issues, it is natural 
to investigate the possibility to extend Whittle-like procedures to a spherical 
wavelet framework, so to exploit the double-localization properties (in real and 
harmonic space) of such constructions. This is the purpose of this paper, where 
we shall focus in particular on spherical needlets. 

Spherical needlets are a form of second-generation wavelets on the sphere, 
introduced in 2006 by 01] and 05], and very extensively exploited both in the 
statistical literature and for astrophysical applications in the last few years. 
Stochastic properties of needlets when used to estimate spherical random fields 
are developed in [5J, [BJ, [7] [3T], [32] and 0T]. Needlets have been generalized 
either to an unbounded support in the frequency domain (Mexican needlets) 
by [10] j |2U an d |22] i an d to the case of spin fiber bundles (spin needlets, 
see [17], and mixed needlets [18]), again in view of Cosmological applications 
such as weak gravitational lensing and the so-called polarization of the Cos- 
mic Microwave Background (CMB) radiation. Concerning the latter, applica- 
tions to CMB temperature and polarization data are presented for instance 
by [5J, [5J, [12], [15], Q]5], QS], 0DJ, 05J, [5D], [52J. Indeed, as described 
for instance in [10] and [5J, satellite missions such WMAP and Planck (see 
htt p://map.gsfc.nasa. gov/[ ) are providing huge datasets on CMB, usually as- 
sumed to be a realization of an isotropic, Gaussian spherical random field. Pa- 
rameter estimation has been considered by many applied papers (see [24] for 
a review), but in our knowledge until now no rigorous asymptotic result has 
so far been produced on these procedures. We refer also to [5J, [IS], [IS], 05] . 
05] . [35] for further theoretical and applied results on angular power spectrum 
estimation, in a purely nonparametric setting, and to [27] . [25] , [25], [3D], [2D] . 
[25] . [34] . [3T] and [35J for further results on statistical inference for spherical 
random fields or wavelets applied to CMB. 

As mentioned earlier, the asymptotic framework we are considering here is 
rather different from usual: we assume we are observing a single realization of 
an isotropic field, the asymptotics being with respect to the higher and higher 
resolution level data becoming available. In view of this, our paper is to some 
extent related to the growing area of fixed-domain asymptotics (see for instance 
[3J, [3BJ); on the other hand, as for [T3J some of the techniques exploited here 
are close to those adopted by [51] . where semiparametric estimates of the long- 
memory parameter for covariance stationary processes are analyzed. In terms 
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of angular power spectrum behaviour, we shall also indeed for semiparametric 
models where only the high-frequency/small-scale behaviour of the random field 
is constrained. In particular, we consider both full-band and narrow-band esti- 
mates, the latter entailing a slower rate of convergence but allowing for unbiased 
estimation under more general circumstances. 

The plan of the paper is as follows: in Section [2] we will recall briefly some 
well-known background material on needlet analysis for spherical isotropic ran- 
dom fields; in Section [3] we will introduce and motivate the Whittle-like mini- 
mum contrast estimators, while in Section [4] we shall establish the asymptotic 
consistency results. Asymptotic Gaussianity is considered in Section [5j while 
in Section [6] we present results on narrow band estimates. Some Monte Carlo 
evidence on performance and comparisons with the procedures in |13j are col- 
lected in Section [7] while some auxiliary technical results are collected in the 
Appendix. 

2 Spherical Random Fields and Angular Power 
Spectrum 

It is a well-known fact in Fourier analysis that the set of spherical harmon- 
ics {Yi m : I > 0, m = represents an orthonormal basis for the space 
L 2 (§ 2 ) , the class of square-integrable functions on the unit sphere (see for in- 
stance [2], [S3], [25], [3S], for more details, and [35], [37] for extensions). Spher- 
ical harmonics are defined as the eigenfunctions of the spherical Laplacian A52, 
e.g. As2Yi m = —1(1 + l)Yi m , see again [S3], [SSJ and [3J] for more discussion 
and analytic expressions. The spherical needlets ([33], [45]) are defined as 

(*) = V^JkJ2 b E (*) Y lm fefc) , (1) 

/ ^ ' m=-l 

where {C^} is a se t of cubature points on the sphere, indexed by j, the resolu- 
tion level index, and k, the cardinality of the point over the fixed resolution level, 
while \jk > is the weight associated to any £ - fe (see also e.g. [7J) and [55]). 
Let Nj denote the number of cubature points for a given level j; as discussed 
by (|44j. [35]). cubature points and weights can be chosen to satisfy 

X jk « B" 2j , Nj w B 2j , (2) 

where by a w b, we mean that there exists c%, C2 > such that c±a < b < 02a. In 
the sequel, for notational simplicity we shall assume that there exists a positive 
constant Cb such that Nj = cbB 2: > , for all scales j. In practice, cubature points 
and weights can be identified with those evaluated by common packages such 
as HealPix (see for instance [5J, [TT], [25]). 

Viewing Lt({x,y)) = Y! m =-i Y i™{x) Y im(y) as a projection operator, 
can be considered a weighted convolution with a weight function b(-) , chosen 
so that the following properties holds (see [33], [45]): for fixed B > I, &(•) 
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has compact support in \B l ,B] and therefore ^(57) has compact support 
in I £ [-B J_1 , ; this implies that needlets have compact support in the 

harmonic domain. Moreover, b(-) £ C°° (0, oo), which will be pivotal to prove 
the following quasi-exponential localization property (see [44J: for any M — 
1,2, ... there exists cm > such that for any x £ E> 2 , 

I , I < 

\rjk — / / ; \\s M ' 

(l + ^arccos((x,^. fc ») 
Finally, we have the so-called partition of unity property: for I > B 



j>0 



which allows to establish the following reconstruction formula (see again |4"4"]): 
for f £ L 2 (§ 2 ) , we have, in the L 2 sense: 

where 

Pjk = (/> ^ifc) L2(S 2) = J ^ ^jk ( x ) I ( x ) dx ■ ( 3 ) 

Consider now a zero-mean, isotropic Gaussian random fields T : § 2 x ft — > K; 
we recall also that for every g £ SO (3) and x £ S 2 , a field T (•) is isotropic if 
and only if 

T{x)=T (gx) , 

where the equality holds in the sense of processes. It is again a standard fact 
(see e.g. [25], [39]) that the following spectral representation holds: 

1 

T ( x ) = ^2 ^2 aimXlm (x) , 
;>o m=-l 

ai m = IT (x) Yi m (x) dx . 

Note that this equality holds in both L 2 (§ 2 x Q, dx (8 P) and L 2 (P) senses for 
every fixed x £ § 2 . For an isotropic Gaussian field, the spherical harmonics 
coefficients aim ar e Gaussian complex random variables such that 

E(azm) = , E(a/ m a iimi ) = S^d^Ci , 

where the angular power spectrum {Ci , I — 1,2,3,...} fully characterizes the 
dependence structure under Gaussianity. Characterizations of the spherical har- 
monics coefficients under Gaussianity and isotropy are discussed for instance by 
[3], [35]; here we simply recall that: 
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Hence, given a realization of the random field, an estimator of the angular power 
spectrum can be defined as: 



1 * 

m— — I 

the so-called empirical angular power spectrum. It is immediately observed that 



E 



m=—l 

Now recall that the needlet coefficients can be written as 

l 



Pjk = 6 ( W ) Y a lrnXl m (£jfc) , (5) 



.B3 , 

Bi- 1 <KBi+ 1 x ' m=-l 



where 



e (p jk ) = ^ k Y, b ( jn ) Y Y ^ fefc) E (o«m) = o • 



SJ'- 1 <i<-B3 + 1 m=— / 

As in [5], we introduce the following regularity condition on the angular power 
spectrum: 

Condition 1 (Regularity) The random field T(x) is Gaussian and isotropic 
with angular power spectrum such that: 

Ci = r ao G(l) > 0, (6) 

where Cq 1 < G(l) < Co, ao > 2, for all I g N ,and for every r G N there exist 
c r > suc/i £/iai: 



< c r u 



for u € (0, oo) 



Condition [T] requires some form of regular variation on the tail behaviour 
of the angular power spectrum C;. For instance, in the CMB framework the 
so-called Sachs- Wolfe power spectrum (i.e. the leading model for fluctuations 
of the primordial gravitational potential) takes the form |6]) , the spectral index 
ao capturing the scale invariance properties of the field itself (a is expected 
to be close to 2 from theoretical considerations, a prediction so far in good 
agreement with observations, see for instance [TU] and [33 J. In particular, this 
Condition will be necessary to prove needlet coefficients ^ to be asymptotically 
uncorrelated (see [5] ) . For asymptotic results below, we shall need to strengthen 
Condition [I] as in |13j . imposing in particular 
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moreover 



Condition 2 Condition^ holds, and 

G(l) = G a (l + 0(r 1 )) , 

whence 



E 



(a,) =c r a ° (i + o^- 1 )) 



As we shall show, Condition [2] is sufficient to establish consistency for esti- 
mator we are going to define. We shall also consider two further assumptions, [3] 
(which implies [2]) , to derive asymptotic Gaussianity, and [4] (which implies [3| to 
provide a centered limiting distribution, see also [13] for related assumptions. 

Condition 3 Condition^ holds and moreover 

G (0 = G (1 + nl- 1 + O (r 2 )) . 

Condition 4 Condition^ holds and moreover 

G(l) = G (l + o(r x )) . 

Example 5 Condition [7] is satisfied for instance by 

0(1) = {'osD'ff, 

where 

P{1) = po + pil + P2I 2 + p m l m ; 

Q (0 = 9o + Qi I + Q2l 2 + q m l m 



en 



are two finite order polynomials such that P (I) ,Q (I) > 0. Condition^ is th 
fulfilled for 6 = 0, k = p m _i/p m — q m -i/q m , while Condition^ holds for n = 
(see also J^,US,JT^). 

Under Condition Q] we have: 

CoB (2-a)j < ^ b 2 ((x,y)) < Cl B^. (7) 

Indeed 

= c [ b 2 (x)g(x)x 1 - a dx + o :j (l) ^ B {2 - a)l ; 
Js 2 

more details can be found the Appendix, Proposition |24| 

As mentioned before, in [5] it is proven, in view of (l7j), that: 
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Lemma 6 Under Condition^ there exists M > such that 
\Cor {P jk ,f} jk 



Cm 

(l + BH{{Z jk ,i jk ,))Y 



< 



As a direct consequence of this lemma, needlets coefficients at any finite 
distance are asymptotically uncorrelated, and hence asymptotically independent 
in the Gaussian case. 

Following ([5]) and [44] . we have easily that: 

Eft = E b2 (^j) E 

k BJ- 1 <l<Bi + 1 V ' m=-l 

= E * (^Wds,, 



Bt- 1 <l<B3+ 1 



where, by equation Q: 

E (E^) = E * (±) (21+1)^(1 + 0(1-1)) . (8) 

The following Lemma provides the asymptotic behaviour of Var (^2 k @j k )- 
Lemma 7 Under Condition^ we have 



lim , 



rW Far |E^J= G ^ 2 («o,5) 



and 



1 r «i ATjj ^ 



= G§T*(a ,S) ,/orjx = j + l 

= G 2 T 2 _(a ,B) , for 31= 3-1 
= ,for lix-il >2 , 



where 



f B 

a 2 (a ,B) : = 4/ b 4 ^)^ 1 " 2 " ^ , 
is- 1 

r 2 (a ,5) : =4/ 6 2 (a;)6 2 (^)x 1 - 2a Vx 



--(«o,S) : =&=*r ] i b\x)b\-)x^ dx. 



4 - s 
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For the sake of notational simplicity, in the sequel we shall write a 2 , t\ . t 2 _ 
(omitting the dependence on ao,B) whenever this does not entail any risk of 
confusion. 

Proof. Simple calculations based on Q lead to: 
and, for ji < j 2 , 



-y2 

J l ' 



fc 1 = l fc 2 = l 



E E fe2 h^r 2 ^ )(2h + l)(2l 2 + l)Cov{d h ,C l2 } 

= J E££U ^(bW)^(bW) (2i + 1) 2C? , for j 2 = j t ± 1 
I , otherwise 

Finally, we have: 

l im 2 y ^f—) (21 + 1) G 2 

l^Bi- 1 V 7 

i \ i r 2ao i 



= 4G 2 , lim V 6 



7— »-oo 



G ° ( 4 / S 64 ( x ) xl " 2Qoda; ) =: g o ct2 ' 



and, if j 2 = j + 1 



= GI\a[\ 2 (x) b 2 (I) x 1 " 2 ^ =: G 2 r 2 



while if j 2 = j — 1 

S3 



lim 2 y fe2 f J_\ & 2 f ' ^ (2| 1} Q 2 

Go f 4 y fo 2 (as) & 2 (Bx) x^^dx^j = 

J\ 2 {x)b 2 (^)x^dx\ =: G 2 r 2 _, 



- G 2 4 
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as claimed. ■ 



3 A Needlet Whittle-like approximation to the 
likelihood function 

Our aim in this Section is to discuss heuristically a needlet Whittle-like approxi- 
mation for the log-likelihood of isotropic spherical Gaussian fields, and to derive 
the corresponding estimator. We start from the assumption that needlet coef- 
ficients can be evaluated exactly, i.e. without observational or numerical error, 
up to resolution level Jj,. This is clearly a simplified picture, analogous to what 
we assumed in [13| for the case of spherical harmonic coefficients; however in 
the wavelet case the assumption can be considered much more realistic. Indeed, 
it is shown for instance in jS] that the effect of masked or unobserved regions 
is asymptotically negligible, in view of the localization properties of the needlet 
transform. Hence we believe our results provide a useful guidance also for re- 
alistic experimental situations. Needless to say, the maximal observed scale Jl 
grows larger and larger when more sophisticated experiments are undertaken: 
indeed Jl is a monotonically increasing function of the maximal observed mul- 
tipole L. The latter is for instance in the order of 500/600 for data collected 
from WMAP and 1500/2000 for those from Planck. In terms of our follow- 
ing discussion, it is harmless to envisage that B Jl+1 = L. The analysis of 
frequency-domain approximate maximum likelihood estimators based on spher- 
ical harmonics is described in [13] , while narrow-band, wavelet-based maximum 
likelihood estimators over R can be found in [43] . 

To motivate heuristically our objective function, consider the vector of coef- 
ficients 

Under the hypothesis of isotropy and Gaussianity for T, we have that j ~ 
iV(0,r), where 

T=[Cov{P ]k ,P 3k ,)] kk , = v (i) • 

In view of Lemma [6] and equation ([2J), it is to some extent natural to consider 
the approximation 




where If?, denotes the Nj x Nj identity matrix. We stress, however, that the 
present argument is merely heuristic - indeed, for instance, elements on the first 
diagonal do not converge to zero. The approximation however motivates the 
introduction of the pseudo-likelihood function: 

C = (2tt)-^ (detr)"' exp f-i^r" 1 /) 
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and the corresponding log-likelihood as: 



E 



log 



NJ 1 Ei b 2 Oj) (21 + 1) C0) {N- 1 & 2 (^) (21 + 1) cm 



up to an additive constant. The full (pseudo-)likelihood is obtained by combin- 
ing together all scales j, so that 



Jl 

E 



E lo § 



1% 



N- 1 Eib 2 ( 1 ^j)(2l + l)Cm 
Let us now introduce the following: 
Definition 8 For a G (2, +oo = A, define the function 



K 



DJ 



with derivatives Kj u (a) given by 



K jtU (a) = (a) 



(-iy 



(2i + i)r a , 



i 



E> 2 ( nj + (bgO u . 



Our objective function will hence be written compactly as: 
Kj L (G,a) : = -21og£(G, a; ^) 



Jl 

E 
E 



(a) 
1 Sfc 



E^ 



G^ (a) 



More precisely, in view of Condition[2]and the discussion in the previous Section, 
the following Definition seems rather natural: 



Definition 9 The Needlet Spherical Whittle estimator for the parameters (ao, Go) 
is provided by 

(aj L i Cj L 



arg min TZj r (G.a) 
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Remark 10 To ensure that the estimator exists, as usual we shall assume 
throughout this paper that the parameter space is a compact subset of M 2 ; more 
precisely we take a € A = [01,02] , 2 < a\ < 02 < 00, and G 6 Y = [7 1; 7 2 ] , 
< 7 X < 7 2 < 00. This is little more than a formal requirement that is standard 
in the literature on (pseudo-)maximum likelihood estimation. 



We can rewrite in a more transparent form the previous estimator follow- 
ing an argument analogous to [5T], i.e. "concentrating out" the parameter G. 
Indeed, the previous minimization problem is equivalent to consider 



(aj L1 Gj L ^ 



argminftj,, (G,a) . 

a,G 



It is readily seen that: 
d 



Jl 



dG 



K Jl (G,a)=J2 

3=1 



1 Z k 0; 



G 2 Kj (a) G 



d_ 
dG 

and because 



K,, L (G, a) = G = G (a) := 



Jl R 2 
L~,k=l Pjk 



E 



(9) 



, Jl N 3 

i£E 



3 = 1 k=l 



2J%_ 
G Kj (a) 



we obtain 
d 2 



dG 2 



Uj L (G,a) 



G=G(a) 



, Jl 

G^fUti 

5£^>o. 

G(af 



G (a) Kj (a) 



Hence G (a) maximizes 1Zj l (G, a) for any given value of a. It remains to 
compute 

(10) 



where 



Rjl (") 



a = argmin {Rj L (a)} 



n 



Jl 



(G(c 



- 1 



log 



3 



Jl 



^ Ka (a) V L N- ^ 
j=i 3 x ' 1^3=1 iv J 3=1 



^NjlogKjia) 
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4 Asymptotic results: Consistency 

In this Section we investigate the asymptotic consistency of the estimators aj L 
and Gj L . In particular we will apply a technique developed by [5] and [5T] , see 
also [13] . 

Theorem 11 Under Condition^ as Jl — > oo we have: 

(uJliGjA -> p (a ,G ) . 

Proof. Let us write: 

AR Jl (a, a Q ) := R,j L (a) - Rj L (a ) 

, G (a) , G fa ) 1 A Ar , If* (a) , G (a) 



where 



£/,/ L (a, a ) - Tjr £ (a, a ) 



G(a) : V-jv ^M 

rp f v , G(a ) G» 

T ^ (a ' ao) 1 =1 ° g GM" l0g G(a) 

r> , , 1 , if,- (a) , G(a) 



It is easy to see that: 



G (oto) = G , log -g^ = log =±- £ N 3 §^ 
G (a ) Ej ^ -Kj («) 



The proof is then completed with the aid of the auxiliary Lemmas 14 15 
which we shall discuss below. In particular 



Pi(\aj L -a \ > e) < Pr inf ARj L {a,a ) < 

\\a— a \>e 

< Pr( inf [Cr Jj .(a J o )-r Ji .(a l ao)]<0 

\Ja— ao|>e 

The previous probability is bounded by, for any 5 > 



Pr inf Uj L {a,a )<S)+Fv\ sup Tj L (a,a )>0 \ ; 

\Jat-a |>e / \|q-q |>£ / 
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for ckq — a < 2, it is sufficient to note that 



lim Pr sup Tj L (a, do) > 



from Lemma 15 
such that 



while from Lemma 14 there exists S £ — g2 + g _ 1 + ^_f 1 log B > 



lim Pr inf C/j^ (a, a ) < <5 e = 

L— >oo \|a — Q |>e 



For ao — a = 2 or ap — a > 2 the same result is obtained by dividing ARj L (a, agj 



by, respectively loglog-B Jl/ or logS" 71 - and then resorting again to Lemmas 
15 Thus aj L — > p ctQ is established. 
Now note that 



14 



G(a jL )-G ( 



i g E/, < 



,2 



1 ^ gpKj (Op) 

Y, j N j j? 1 K j (aj L ) E. ; ,\ , £j ^(ao) 



By adding and subtracting / (£?, ao, <5j L ) ~~ L where / (B, cxq, aj L ) is defined as 
( 28 ) in Proposition 24 we obtain: 



G(aj L )-G 



< 



Kj (a ) ( E fc /3 



1 JL 



+ 



+ 



Jz. 



(«o) 



I (B,a ,aj L ) 



— i— ^Tg (I(B,a ,aj L )-l) 
iy i ,-=i 



|G A | + |G B | + |G C | ■ 



By Proposition 24 we have that 



Kip) =B i{^- ao ) I{B ^ )+ {X) 
Kj (aj L ) 



Clearly 
Pr{|G,|>|} 



< Pr 



\Ga\ > 



n 



\ao - aj L \ < 



Pr <^ \a - ctj L \ > 



< Pr 



Go 



E,- Nj 



ato, otj L 



) B j ( Sj - 



ao) 


Efe P]k 1 


> e 


) 




G Kj (ao) 







oj L (l) 
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< 



< 



1 Go 



aj L- a °)E 



<3 3 j 
1 Go 



aj L -ot ) 



fWj G Kj (a ) 



rm GqKj (a ) 



o.J L (1) 

1/2 



G Go J (-B, ap, SjJ /( a - ao ) +1 )j L 
e B^ 



oj L (l) 



e 

in view of the consistency of a j t . As far as Gb is concerned, we obtain, for a 
sufficiently small 5 > 0: 



Pr{|G fl |>|} = Pr{[|GB|>|]n[logB^> -a 7i |< 



Pr 



= Pr 



{[|G 
{[\Gb\ 



> 



> 



n [logB 1 \a -aj L \ > 
n [log B J |a Q - 3jJ < S] } + 0ji (1) 



Because for < x < 1, we have |e 21 — 1| < x, we have: 



= |exp(-j (a - aj L ) log B) - 1| 
< (a-o - aJz,)logS J , 



and hence, in view of 

Pr{ |G 



> 



< Pr 



Jl 



3 

Kj (a ) 



n [log^'lao-a^l < 



e 

> - 
~ 3 



< 



C I (B,a ,aj L ) 



3 1 3 = 1 



J2G \ogBm\a -aj L \ 



n [logS J |a -aj- t | < 5] 
C SJ L 



Finally, in view of (29), 



1 Ji 
— G ° ( J ( B > "o, uj l 

1^3 ^3 4=1 



)-l) 



< 



JlGj |a - SjJ 



Hence, because for a sufficiently small 5 > 0: 

Pr{[|G c |>|]} < Pr {[\G C \ > ^] n [\a -aj L \< 6]} +Pr{\a -& JL \> 5} 
= Pr { [\G C \ > J] n [N - 3jJ < <5]} + o 7i (1) 
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Pr { [ |Gcl - II n [K - 3ji 1 < ^} - J^W- = 0J - (1) 



as claimed. ■ 

Here we present the auxiliary results we shall need on G and its derivatives. 
We introduce 



Go (a) 
Gr (a) 
G 2 (a) 



X^ Jl m, ^ Ka (a) 



f/ ,j (a) ; 



-G(a) 



'l A t 



Ui,j (a) ; (11) 



da 2 



where: 



Also, let: 



, U 2;j (a) = 2 



K j> x(a)\ Kj, 2 (a) 



K 3 (a) J Kj (a) 



Go (a) 
Gi (a) 
G 2 (a) 



ri i \ 1 at GoK 3 (a ) 

G (a) = —j—r ^ Nj t/ 0j (a) ; 



(/a 



V Jl - AT- ^ Ka (a) 
Gin) ' v >:.v/ / '; A ^" : r l ,12) 



d 2 
da 2 



Jl 



nl \ 1 i\r GqKj (ao) rr , \ 



The first result concerns the behaviour of expected value and variance of the 
estimator G (ao) computed in ao, the second regards the uniform convergence 
in probability of the ratio between G and G, and their fc-th order derivatives 
Gfc, Gfc. 

Lemma 12 Let G (a) be as in |$|). Under Condition^ we have 



E 



(G(a Q ) 



Go 



limB 2jL Var(G(a )) = Glp 2 (a ,B) 



B 2 -l 
B 2 ' 



where Iq (B) is defined by (21) in Proposition 24 and 



p 2 (a ;B) 



a 2 (a ;B) + B- a °T 2 (a ;B) 
I 2 o(B) 



t 2 (a ; B) := t\ (a ; B)+t 2 _ (a ; B) 
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Proof. By ((8), we obtain that 



E 



(6(a)) 



AT ^-^ 



= Go + 0,7,(1) 



^7 (ao) 



while from Lemma [7] and Proposition 24 (see also the proof of Lemma 18 1, we 
have 



Var (5(ao)) =G^ 2 (a ,S) 



g 2 -l 
B 2 



B~ 2Jl + o(B~ 



2 Jr. 



as claimed. ■ 

Lemma 13 Under Condition^ we have for n = 0, 1, 2: 

G„ (a) 

sup r i \ ~ 
(*n (a) 

Proof. Under Condition [3j we observe that: 



G n (a) 
G„ (a) 



1 = 



V Jl - N- GoKj( - ao ' ) U- (n) 



E Jl fW K i\ a o) TT { 
,=! V JV J ir,( a ) U 3,ri \ 



Then we have: 

2jj=l V iv J KTfal ^ a 



^"Efc (g 



< 



V 



V ,/W- Ki{aa] U- (a) 



sup 



('■to) 



7 2 



>5 F 



>6 F _ 



J 



In view of Proposition [24| and equations ( 34 ) and ( 35 ) in Corollary 25 de- 
scribed in the Appendix, we have 



U Jt i (a) 



U j<2 (a) = 2 



Kj (a) 



\ogB J +Oj{l) 



{K. hl (a)) 2 K ji2 (a) 
{Kj («)) 2 Kj («) , 



= (log^) +o i3 (l). 
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Then. 



o{b- Jl ), 



so that 



sup 

L 



< +00 



E&^§^>(«) 

Also, by Markov inequality and Lemma [7], we have that, for all j: 



/AT. 



0- 



jk 



- 1 



>CJ 2 L 



< 



CJl 



( 1 @jk 

1 l -Var(y ^ 

CJlGlN 3 Var 



- 1 



CJl 



Kj (a ) 



whence 



sup 

j=i,...,J L 



< Jl sup 

3=1,. ...Jt 



^E 



P]k 



'Nj ^ G K 3 (a ) 



^E 



Go-Kj (ao) 
< J L xO(j- 2 ) =0{J1 1 )=o Jl {1) . 

Under Condition [2] we have: 



p 2 (a ,B) = O(Jz'- 



> CJ L j 
- L >CJ L 



Gn (a) 
G„ (a) 



V Jl - , /W Kj< - ao ^ U- (a) 

2^.7=1 V JV J if j (a) U J.« V«J 



'j\r. ^ 



f l 2 . e[e/&] 



Z-fj=l iv J jr., (a) V a 7 



From ( 36 ) , it is easy to see that the second term in the last equation is O (£> Jl ) , 
while for the first term we follow the same procedure already described, also 



considering that, by ( |36 1 

^E 



ftjk 



E [0 



jk\ 



"J fc 



G Kj (a ) G Kj (a Q ) 



E[0 



fjTj G Kj (a ) \E 
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<a£(i + o(b-'-)) 



We can establish now the asymptotic behaviour of Uj L (a, cxq), for which we 
have the following 

Lemma 14 For all e < — a < 2 

lim U Jl (a,a ) 

Jl—^oo 



Jl 



log 



5 2 -l B 2 (a-a ) 

B 2+{a-a ) _ I B 2 -1 

Moreover, if a — a>o < —2, we have 



logB >5 e > 



lim 

Jl^oo log _B Jl - 



and if a — ao — 2 



Kj (a) 



cvq — a 



1 > 



lim 

J L -yoo log 



(a ) 



Proof. By recalling = cbB 2j (see ([2])) and (26 1, we observe that, 



B?JlB 2 



3=1 



+ log (J(B,a ,a)) + o Ji (l) 



B 2 - 1 
B 2J ^B 2 



Jl 



{a -a^Y^B 2 ' log Bi +o Jl {\) 



Z> 1 r (a-a )logs(s 2Ji L 



B 2J ^B 2 



B 2 - 1 



1 



(a - a ) log 5 



B 2 - 1 



B 2 - 1 
+ o Jl {\) , 



(13) 
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Now, for a — ao > —2: 



using (38) in Proposition 27 described in the Appendix, with Ji = 1 and s = 2. 

Kj (ota) 



log 



3 K s (a) 



log 



log 



1 J-^ .^(2-ao)j 



^.7 = 1 .7 = 1 



+ log (J(B,a ,a)) + 0^(1) 



_ ^ g(2+a-a )J L j^2+a-a 



Q2J L g2 Q2+a-u a _ I 

r b 2 -i 



+ log (J (B,a ,a)) + o Jt (l) 



-ao) — 1 ' 



B ( a - ao )JL Ba - a « [ + iog(I(B,a ,a)) + oj t (l) 



( B 2 — 1 1 

lo § { B 2+(a-a ) _ i j + ( a - «o) { Jl + 1} log B + log (/ (B , a , a)) + oj L (1) . 

(14) 



Hence, combining ( 14 ) and ( 13 ) we obtain 



log 



B 2 - 1 



g2+(a-a ) _ I 



+(a- a ){J L + l}logB-(a-a ) \ Jl 



B 2 - 1 



logB+oj L (l) 



B 2 - 1 1 N (a - a ) 

li >■" < B 2 +[a - ao) 1 ) + (« " «o) log B + log B + 0ji (1) 

b 2 — i i r b 2 i 

= log 4 „,+r„_,^ J + (»-»o)| |> log B + o./ L (1) • 



£2+(a-a ) _ I 

Now consider the function 



{B 2 — 1 1 f B 2 1 

it is readily seen that for x > —2, l(x) is a continuous function such that 

„, x B 2+x \ogB B 2 , 

Or) = ^ — + ^ logB , 

v ; B 2 +z — 1 B 2 — 1 6 ' 

i'(0) = , l'{x) < for x < , Z'(x) > for x > , 

whence Z(0) = : is the unique minimum, and l{x) > for all x ^= 0. The first 
part of the proof is hence concluded. 
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Take now a — ao < —2; we have 
J 



1 



log B 2J t 



log B 2J t 



loe 



log 



Kj (a ) 



1 % 



Kj (op) 



3=1 



log B 2 ^ + Jfi (1) - £ A/j log fl* 



X^ JL AT- f 



ao — a 



1 + 0.7,(1) 



Finally, for ao — a = 2 we obtain 

J. 



1 



log J L 



log 



gj (go) 



gj («o) 



log J 



whence 



J— J - log 5 2 ^ + log J L + O jL (1) - 1 ]T ATj [log + Jl (1)] 1 



— !— {-log B 2J - + log J z + Jl (1) + log S 2J - + Oj L (1) } , 
logJ L 



if, 1 »r K i («o) 1 »r , ^7 (ao) 1 



lim 

Jl^oo log 



as claimed. ■ 

Now we look at Tj £ (a, ao). From (12), we can prove the following: 
Lemma 15 As Jl — > oo, we have 

sup\Tj L (a,a )\ = o p (l) . 



Proof. From (11) and (12), we have that 



G(a ) = 

G(a ) J2jLi Nj G oK 3 (a ) 
From Lemma 1121 it is immediate to sec that 



E 



( G(a ) 
\G{ao) 



1 =0 
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and 

By applying Chebichev's inequality, we have 

G(« ) 



0(B- 2J -) 



G(a ) 



- 1 



and from Slutzky's lemma: 



log 



G(a ) 



On the other hand, in view of Lemma [13 

G(a) 



sup 



G(a) 



- 1 



so the proof is complete. 



, 



> P 



5 Asymptotic properties: The Central Limit The- 
orem 

For the Central Limit Theorem, we shall exploit some classical argument on 
asymptotic Gaussianity for extremum estimates, as recalled for instance by |46j , 
see also [13]. We shall in fact establish the following 

Theorem 16 Let aj L = arguiin^^ R^(a). 
a) Under Condition^we have: 

B Jl 

(a jL - oo) - O p (1) ; (15) 



\ogB- ,L 

b) Under Condition^we have 
B Jl 
log B Jl 



(aj L - ato) — > P 



where 



B + l 

m = - k ; 



whe 



z) Under Condition^we have: 

B Jl -(aj L -a )^Af(0,D) 

re 

D = D (a , B) = p 2 (a ; B) (B) , * (B) 



(16) 



B 4 log 2 B 



(17) 
(18) 
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Remark 17 In particular, the previous result implies 

(a.j L - a ) = O p (^y^) , as L -t oo , 

showing that the rate of convergence of Whittle needlet estimates is the same as 
for the spherical harmonic estimator analyzed in ilStf . 

Proof. By a standard Mean Value Theorem argument and consistency, for 
each Jl there exists aj L € (ao — a.ao + a) such that, with probability one: 

/- s Sj L (a ) 

(a jL - a ) = - — — - , 

Qj L ( a L) 

where Sj L (a) is the score function corresponding to Rj L (a), given by: 

1 W ^ k ( K i+( a ) \ I 1 ^ %( a ) » r 

Ej ^'YV^ (») K i («) v A >" > 3 

1 G(a) v / %(«) \ V /G(«) G(a)^ 



Ej G (a) Y V K i ^ ' k V G («) G (a) ^ (a) G (a) , 



and 



i.e. 



d 2 



G 2 (a) G (a) - {G 1 (a)) 2 1 ^ K ji2 (a) K 3 (a) - (K jA (a)) 2 
(G(a)) LjMj . (Kj(a)) 



at Kj{a a ) I ( Kj A (a) \ 2 K ji2 (a) \ \ ( tr- AT K j( a o) 



Ki(a) J) t 1 ^ Ar # ii2 (a) (a) - (K jtl (a)) 



3 



where G(a), Gi(a), G^a) are respectively the estimate of G and its first and 
second derivatives, as in Lemma [l3j In order to establish the Central Limit 
Theorem, we analyze the fourth order cumulants, observing that this statistics 
belong to the second order Wiener chaos with respect to a Gaussian white noise 
random measure (see [47] )■ Let 

B J -Sj L (a ) = ^ JZ J2( A i + B ^ > 

3 
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where 



B: 



B 2j \ogB j 
B 2j logB J 



NjGoKj (a ) 

Gj L {aq) _ 
Go 



- 1 



(19) 
(20) 



In the Appendix, Lemma 29 shows that: 



^™ | + %), J2(A h + B j2 ), J2(A j3 + B j3 ),J2( A M 

\ li h h 1-4 



o 



Jl 



J 4 L log 4 B 

B 2J L 



Central Limit Theorem follows therefore from results in [47 . The proofs 
of (16) and (17 1 are completed by combining the following Lemmas 18 and 
19 Observe that under Condition [3j the only difference between (15) and (16) 



concerns the possibility to estimate analytically the bias term. 
The following result concerns the behaviour of Sj L (a). 



Lemma 18 Under Condition^ we he 

B j l 
\ogB J i ' 

while under Condition^ we have: 



B + l 
B 



-k ; 



B J ^Sj L (a )^Af(o,p 2 B (a 



log 2 B 



Proof. We have that: 



Go 



G(a) 



K j ( a o) 7 \G (a ) G Kj (a ) G 



where we recall that for Lemma fl3] 

Go 



G(a ) 

Then we will study the behaviour of 



s Jl (« ) = 



Kj (a ) 



E 



G Kj (ao) 



G(a ) 
Go 
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Under Condition [3j simple calculations, in view of (37), (34) in Corollary 25 
and (39) in Proposition 27 lead to 



lim 

Jl— too 



r log b Jl 



= lim 

Jl— >-oo 



lim 

Jl— >oo 

B + l 



log B ,,L 
B J t 



E(S Jl (a )) 



_^L_£i?;iogij'+ 0j£ (i) 



while under Condition 3] we have 



lim E (Sj L (a )) = 

Jl->°o 



Moreover we obtain: 



Var (S Jl (a )) 



where 



B 



C = 



Var 



1 >p g^i (go) / Sfc/^jfc _ „ G(ao) \ 



j 

= A+B+C, 



EjNjf^ K^ao) \GoKjiao) 3 G J 



A = 1 \^ ( Kj U l (g0) K j 2 ,l ( a 0) \ C ( Efcl Pjlfcll Sfc 2 1^72*2 



Ei N j fa 



K h (a ) K j2 (a ) / \ Go-ftT^ (a ) ' G i^ 2 (a ) 



_J y- f ^M^W ) iV . ^ Var ( gw) . 



-2 
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In view of Proposition 24 and Lemma [7J we obtain 
A - 1 y ( K ^ l{ - aa) \ 2 Var ( ^ k ^ k 



Jl-i 

E 



(«o) 



(ao) («o) 
Kj (a ) K j+1 (a ) 



Cow 



Efc ^j* Efc /3. 



G Kj (a ) ' G K ]+1 (a ) 



/ g^i (go) gj-i (ao) \ „ / Efcffjfc Efcfij-i,fc \ 
£j V ^ (oo) tfj-i («o) / ^ G JT, (ao) ' G^-i («o) J 



+ 0J £ (1) 



E 



P 2 (ao,S) 



^ B^' log 2 B-' + o Jl (1) 



because 



- E 

■ E( 

and, likewise, 



A. ^ ^ (ao) Kj+1 {ao) j y 



G iTj (ao) ' G Kj +1 (a ) 
^ (a ) (a ) J / 2 (a ,B) B<*° B~^oU+D 



Oj (1) 



^■,i(a )\ 2 B^r 2 , 



Kjiao) ) I$(aa,B) 



oj (1) 



( X J -i(a )g J -i(a ) \ / _J> 



Efe /3 2 fe E 



Kj (a ) G Kj-i (a ) 



E 



JT,-(a ) 7 / 2 K,S) 



■ 0j -(l) 



On the other hand, by Lemma [l2j we have: 



Jl 



Jl 



Jl 



D 



(E, 



lP 2 (a , B) ^ B 2jl logfl* B * h ^ogB 32 B2J3 + °Jl 



ii=i 



J2 = l 



J3=l 



= P 2 K,S) 



(e,^) 3 



oj, (1) ■ 
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Finally we have: 

C 



(E,iv,)- 



ii j2 



gjij (gg) gja.1 (gg) 

i^-j (a ) K h (a ) 
(«o) ^2,1 ( a o) 



J3 



I G (a ) iiTji (ao) ' G K J3 (a ) 



^ \ K h (ao) if ia (ao) J J2 i 2 (B) 



-2p 2 (a ,£) 



o- 2 + B- ao T 2 ) 



B^+o jL (l) 



o.J L (1) ■ 



Hence, following Corollary 28 in the Appendix and equation (38), we have: 

1 



Var(Sj L (a )) = P 2 (a ,B) 



(E^) 

^ // J ' ^ B 2j log 2 //•' I £ // ' ' log B j I + o (// ' 



= p 2 (ao,B) 



:Z Ji (2)+ (i?- 2 ^) 



s (5 2 -ir 47 „ B 6 
7 ~B AjL \og 2 B— —+o(B 



P 2 (a ,B) 



B eJ ^B e 

.2 



(B 2 - iy 



Finally we have: 



P 2 («o,B)0-^B-^ +o(B-^) 



log 2 B 



(21) 



^ ( B ' ,LSj t ("")) = ^ ( a 0> B ) / S 2 _ ' 

as claimed. ■ 

The following Lemma regards instead the behaviour of Qj L (a). 
Lemma 19 Under Condition^ we have: 

BHog 2 B 

Proof. By using [T3j we obtain: 



(a) = G 2 (a)G(a)-(G 1 (a)) 2 | ^1 y. ^ ^ 2 (a) ^ (a) - (if iiX (a)) 2 



(G(a)Y 



Ei^V 



(a))' 



2G 



L,j^3 Kj ( a ) K 3 { a ) ) K 3 ( a ) j j \2^j^3 K 3 ( a ) ) 

(E^-f$) 2 

,r K {cc ) ( K jA (a)\\ 2 

ly j Kj (a) { K 3 { a ) )) ^ 1 K j>2 {a) Kj (a) - (K jtl (a)) 



^7 Kila) / J 



K 3 {oc) 

Qj L (a) can be rewritten as the sum of three terms: 

Q.J L (a) = Qi (a) + Q 2 (a) + Q 3 (a) 

where: 

QT m (<*) 

«1 W = nden 



Qf en (a) 

^ N K i (ao) / %i(a) \ 2 \ / v wSHUfr /V M2£) L^tlMlV 
2-j iV J /Y 3 (a) ^ ^(a) j y ^j JV J Kj (a) ) JV J /Y 3 (a) ^ K 3 ( a ) )) 

2^7 JV J iY,(a) J 



Q^ um (a) 
" (a) 

(E^i) (%#) 2 ) - (E^§^) foiV; (%C 2 



•32(a)- 



(e^^)e 

Qr n («) 
Q3 (a) = op^y 



(e^SS*) 




(e^-m) 




\2^j^3 Kj ( a ) Kj(a) ) 






E^ 



The next step consists in showing that: 

Q2 {a) + Q 3 (a) = Oj L {l) . 



Using Corollary 25 Q^ um (a) can be written as: 

Q 2 {a) 



(5) 



2 



°J*(1) 
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while Qg 1 "™ (a) becomes: 
Kj (oo) 



g 3 (a) 



(B) 
h(B) 



(B) 



oj L (l) 



so that: 



Ql um [a] + Ql um (a) 
Qf n (a) 



It remains to study Q^ n (a) ; by using (24) and (27 1, we have: 



lim , , 

j^co B 2(2+^)j t Kj{a) 



c 2 B I (B,a ,a) 



= lim 



c|/(B,a ,a) 



2 1+- 



S 2 - 1 



> 



Finally, we prove that Qi —} p fg|~n^ • Using Corollary 
the numerator Q™ um (a) as: 
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we write 



)™" m (a) 
g j ( Q o) at) ( V" N K i ^ 



log S J 



M5) 

MS) 



ifj (a ) 



. log B J 



MB) 
/o(B) 



_ | K j (go) ^ g J ( a °) 



Let s = 2 (1 + ^2^); by applying (42) we have: 



1 



lim — ; — — i / 1 



(a) = lim 



c 2 B I (B,a ,a) 



° K j (") 



Zj l (s) 



= log'B 



B 3s 



(B> - 1) 



.<? B I (B,a ,aY 
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It remains to study Qf en (a) ; by using again (26) and (27): 



1 



-\de7i 



(a) = 



lim 



c 2 B I (B,a ,a) 



c B I (B,a ,a) 



B s 



B s - 1 



Hence 



lim Q (a) 



B 2 ( 1+ ?)log 2 5 



2 ' 



B 2(l+^)_ 1 

For the consistency of Sl, for g [oq — 5l, «o + Sl], we have 

B 2 log 2 B 



Q(a L ) 



P (B 2 - l) 2 



lim Var (B j -S Jl (<*„)) = p 2 (a ,i?) /"f^ . 

■ 

To investigate the efficiency of needlet estimates, fix B Jl = L/B, so that 
the highest frequency covers the multipoles I = [L/B 2 ~\ + 1, L; observe that 

B Jl+1 (aj L -ao)=L (aj L - a ) A N (0, B 2 x D (B, a )) , 
while parametric estimates based upon standard Fourier analysis (see [13]) yield 

L(a L -a Q )^Af (0,8) . 

For any given value of B, the asymptotic variance D (B,ao) can be evaluated 
numerically by means of ( 18 1 and a plug-in method, where ao is replaced by 
its consistent estimate aj L . In practice, though, this is not really needed for 
the values of B which are commonly in use, i.e. B ~ 1.1/1.5. In fact, using 
logx ~ x — 1 + o(l) as x — > 1 we have 



B-vl 5-1 



5 2 



1 



= lim 



1 (B + l) d (B - 1) 



B-H B - 1 



B 2 
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A standard choice for the function b(.) (see [5], [32]) is provided by 
b 2 (x) = b 2 (x;B) = , for x £ (=,B) , 



b\x) 



b 2 (x) 



B 

r(^~ B-i ( x ~ S 1 )) / 1 \j 

-i i-i exp(- CT )rfu i 

1 1 j for — < x < 1 , 

J_ 1 cxp(- TZ U,)du # 

j , tor 1 < x < B . 

J_ 1 exp(- J ^ I )du 



For this choice of &(.), analytical and numerical approximations allow to show 
that 

2 



whence 



lim (B - l)p z (a; B) = 1 , 

B-»l 



Jm D (B, Q0) . ,„ ( < B+ '» + fl f < B + " (1)) - 8 . 

Summing up, the variance of the needlet likelihood estimator is very close to the 
"optimal" value (e.g. 8) which was found by [13] for the Fourier-based method. 
Some numerical results to validate this claim are provided in Table 1 for a range 
of values of B and ao. These numerical results are confirmed with a remarkable 
accuracy by the Monte Carlo evidence reported in Section [7] below. 



B 


^2 


^2 


^2 


2 


a 


2 


3 


4 


2 


3 


4 


2 


3 


4 


2 


3 


4 


o 2 


0.27 


0.27 


0.27 


0.53 


054 


0.54 


1.15 


1.16 


1.16 


2.09 


2.10 


2.10 




0.04 


0.05 


0.05 


0.13 


0.13 


0.13 


0.44 


0.44 


0.44 


0.58 


0.58 


0.58 


Io 


0.17 


0.17 


0.17 


0.35 


0.35 


0.35 


0.70 


0.70 


0.70 


1.39 


1.39 


1.39 




8.46 


8.48 


8.50 


5.00 


5.04 


5.09 


2.58 


2.61 


2.63 


1.40 


1.41 


1.43 




0.75 


1.18 


2.08 


3.51 


B' Z D 


7.67 


7.69 


7.70 


8.36 


8.43 


8.52 


10.7 


10.8 


10.9 


20.6 


20.8 


20.9 



Table 1: Some deterministic results for different values of B and ao- 



Remark 20 It is shown in |5]/ how needlet coefficients are asymptotically unaf- 
fected by the presence of masked or unobserved regions, provide they are centred 
outside the mask. It is then possible to argue that the asymptotic results pre- 
sented here remain unaltered in case of a partially observed sphere, up to a 
normalization factor representing the so-called sky fraction, i.e. the effective 
number of available observations. This is a major advantage when compare to 
standard Fourier analysis techniques - in the latter case, asymptotic theory can 
no longer be entertained in the case of partial observations. Again, for brevity's 
sake we do not develop a formal argument here; proofs, however, can be routinely 
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performed starting from the inequality 

E{P jk -(3* k } 2 < Cm 

W% ~ {l + BH(i ]kl G)} M ' 

valid for every integer M > 0, some constant Cm > 0, for G denoting the 
unobserved region, see again J3f,Jj5f. 



6 Narrow band estimates 

As discussed in the previous Section, under Condition [3j asymptotic inference 
is made impossible by the presence of the nuisance parameter m. It is possible 
to get rid of this parameter, however, by considering narrow-band estimates 
focussing only on the higher tail of the power spectrum. The details are similar 
to the approach pursued in analogous circumstances in [13] . We start from the 
following 

Definition 21 The Narrow-Band Needlet Whittle estimator for the parameters 
d = (a,G) is provided by 



Jl 

(aj L . Jl ,Gj L . Jl ) := argmin ^ 



a,G 
' 3 = Jl 



GKi(a) ^° g lG^(a) 



or equivalently: 

&Jl;Ji = argminiZj 1 .j 1 (o! ! G(a)), (22) 

a 

where J\ < Jl is chosen such that B Jl — B Jl — > oo and 

B* = B* (1 - a {Jl)) ,Jr=J L + . (23) 

log ±> 

We choose < g (Jl) < 1 s.t. lim^-^ g (Jl) = and Iiruj^oo J\g (Jl) = 



For notational simplicity B Jl is defined as an integer (if this isn't the case, 
modified arguments taking integer parts are completely trivial). For definitc- 
ness, we can take for instance g (Jl) = 3 • 

Theorem 22 Let aj L .j 1 defined as in 



'22). Then under Conditionlftwe have 
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where 



4> (B) = log 2 B 



B 2 



(B 2 — 1) V(^ 2 -l) 



log B - 1 



Proof. Because the proof is very similar to the full band case, we put in 
evidence here just the main differences. Consider: 



Qjl;Ji ( a ) 



da 
da 2 



Rjl-Ji (") 



RJl;Ji («) 



Let 
Sj l -j x (a ) 



Sj l -,j 1 (a ) 
-1 



Go 



*- , j=Ji ■> j=Ji 



G(a ) 



G(a ) \ 



N 4 hi K j( a o) fr[\G(a Q )K : j{a ) G (a ) J 



We have: 



E {Sj L . tJl (ot )) = KlogS 



5 (i-^; 



o (J L B- 



B + l 1 B(B-l) , , 
k log 5 — ^— J L ^ / +o(j L g-^ 



B B^ B 2 — 1 



1 



K\ogBj^-J L +o(.J L B 



Following (41| and (21), we have 



VariSj,,^ (a )) =p 2 (a ,B) / (2) 



For (23), we have: 



(2) 
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1 - - (1 - 9 [J L )f (1 - log, (1 - 9 (W)) 2 

^\ 2 

(1-2 5 (J L )) 
x (1 - log B (1 - g (J L ))f +0(g (J L ) 2 ) 

AZ., L , h (g(J L ))+o(g(J L ) 2 ) 



B 2 - 1 + 2g (J L )V _ (B*-l 
B 2 { B 2 



B 2 - 1 



B 2 



B 2 -1V/ 4 / 2 

^B 2 ^) \{B 2 -l) + y~\ogB 



^g(J L )+0(B iJ -g{J L ) 2 ) ,(24) 



where 



1 + lW^) l '' M )-( 1 + {iiB-^ siJL) 
1 



\B 2 -l) V lo S B . 
Thus we have 

Z., L ,, h (2) = B iJ - $ (S) 5 ( J L ) +0 (B iJ -g (J L f 



Note that $ (S) > for B > 1. 

On the other hand, simple calculations on Proposition 27 lead to 

3 



E s2j ) = fB2 g6 n3 (b 2Jl -g 2(Ji ' i) ) 3 + (i? 6j ") 



(B 2 - 1) J 

J0^^ L (l - S- 2 (1 - 3 ( J L )) 2 ) 3 + a (B^) 



B ( > 



^(^ 1 ) 3 + 0(B^g(J L )) 



(B 2 -l) 3 ' 
= B 6Jl +0{B 6J ^g(J L )) , 

hence we have 

Var (Sj l;Ji (ao)) = p 2 (a , B)$ (B) «? ( J L ) £T 2Ji 
Consider now Qj l; j 1 (ck), which we rewrite as 

Qjl Ji (<*) = 7^ • 
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Following a procedure similar to Lemma 19 we have 

Qnum («) = c%Gp (B, a , a) Z Jl . Ji (s) 



where s = 2 (l + SL ^ M )- Following (24 1, we obtain 



Qnum (a) = c 2 B G 2 I (B, a , a) $ (B, s) B 2sJ ^g{J L ) + O [B 2s ^ g (J L ) 2 j , 
where 



$(B, S ) = log 2 i?_ 



2 u_ B L_ ( 2s 9(Jl) slogB-2 
(B s - I) 2 \ B s -I logB 



Finally, we obtain 

Qden (a) = c 2 B G 2 I (B, a , a) I B ° j 

\j=Ji 

= c 2 B G 2 Q I (B, a , a) B 2sJl + o (B 2sJl ) 

Hence 

Qj L ,. h (a) = $ (B, s) g (J L ) + (B 2sJ ^g{J L y 
and, for the consistency of a, we have 

Qj l ;Ji(<x)^p$(B)9(Jl) ■ 

Thus 

/ Qj L ;Ji ("J 

as claimed. Finally we can see that 

g{J L )-B ^^ ^ =0(j L -g(J L yA -> 0. 



Remark 23 A careful inspection of the proof reveals that the asymptotic result 
could be alternatively presented as 



B 6 



B~ Jl y/Zj LiJl (2) (a L , Ll{5) - a ) A AT \0, p 2 (a , B) 
where 

z.,^ (2) = f £ bA ( f; s^ 2 iog 2 s ] - ( f; B^iiogB 
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Fourier estimates 



Needlet estimates 




-4 -2 2 -4 -2 2 



Figure 1: Distribution of normalized (oll — ao) and (3j L — ao), L = 2048, 
a = 2. 

7 Numerical Results 

In this Section we provide some numerical evidence to support the asymptotic 
results discussed earlier. More precisely, using the statistical software R, for 
given fixed values of L, a , B and G and the alternative conditions discussed 
in the previous Section, we sample random values for the angular power spectra 
Ci, and evaluate the corresponding needlet coefficients f3j k ; we implement stan- 
dard and narrow-band estimates with both standard Fourier (as described in 
p"3] ) and needlet methods. We start by analyzing the simplest model, i.e. the 
one corresponding to Condition [4] Here we fixed Gq — 2. In Figure I, the first 
column reports the distribution of Fourier estimates of ul — ao normalized by 
a factor ^/2L/A, while the second column reports the distribution of a.j L — ao 
normalized by the factor D(ao, B)~^B Jl . In Table 2, we report the sample 
means and variances for different values of L and ao , while in Table 3 we report 
the corresponding Shapiro- Wilk test of Gaussianity results. Figure 2 describes 
graphically the behavior of normalized distributions of estimates of ao in both 
classical Fourier and needlet analysis, full band and narrow band, with k = 1, 
under Condition [3] Table 4 provides sample means and variances for different 
values of ao , with = V2, L = 1024, L\ — 724, and the results of the corre- 
sponding Shapiro- Wilk test of Gaussianity. Overall, we believe this numerical 
evidence to be very encouraging; in particular, we stress how the asymptotic ex- 
pression reported earlier provide extremely good approximations for the Monte 
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Carlo estimates of the standard deviation. 



B = 2 






L 


a 


mean 


sd 


(L ■ sd) 2 /B 


mean 


sd 


(L ■ sd) 2 /B 


256 


2 


1.9984 


1.12 • 10" 2 


1.03 


1.9981 


1.68 • 10" 2 


1.12 


3 


2.9994 


1.13 • 10~ 2 


1.05 


2.9997 


1.84 • 10- 2 


1.07 


4 


4.0009 


1.10 • 10~ 2 


0.99 


3.9996 


1.89 • 10-' Z 


1.12 


512 


2 


2.0005 


5.79 • 10" 3 


1.09 


2.0002 


8.55 • 10" 3 


0.93 


3 


2.9995 


5.76 • 10~ 3 


1.08 


2.9999 


8.50 • 10~ 3 


0.91 


4 


3.9997 


5.59 • 10" 3 


1.02 


3.9997 


9.35 • 10" 3 


1.09 


1024 


2 


2.0002 


2.79 • 10" 3 


1.02 


2.0002 


4.42 • 10" 3 


0.99 


3 


2.9999 


3.01 • 10" 3 


1.18 


2.9998 


4.40 • 10" 3 


0.98 


4 


3.9997 


2.82 • 10~ 3 


1.04 


3.9998 


4.39 • 10~ 3 


0.97 



Table 2: Sample means and variances of Ox and aj L , for different values of L 

and ao, B = 2. 





oil - Shapiro test 


aj L - Shapiro test 


L 


a 


W 


p- value 


W 


p- value 


256 


2 


0.9921 


0.35 


0.9952 


0.32 


3 


0.9945 


0.68 


0.9939 


0.59 


4 


0.9906 


0.22 


0.9909 


0.24 


512 


2 


0.9956 


0.83 


0.9945 


0.12 


3 


0.9915 


0.30 


0.9943 


0.65 


4 


0.9946 


0.70 


0.9895 


0.15 


1024 


2 


0.9885 


0.11 


0.9955 


0.82 


3 


0.9967 


0.94 


0.9949 


0.74 


4 


0.9915 


0.29 


0.9959 


0.87 



Table 3: Shapiro- Wilk Gaussianity test of and aj L , for different values of L 

and ao, B = 2. 





a L 


Shapiro test 




Shapiro test 




mean 


sd 


W 


p-val 


mean 


sd 


W 


p-val 


2 


2.004 


2.68 • 10" 3 


0.9985 


0.35 


2.007 


2.75 • 10" 3 


0.9989 


0.32 


3 


3.004 


2.76 • 10" 3 


0.9988 


0.91 


3.004 


2.79 • 10" 3 


0.9988 


0.94 


4 


4.004 


2.88 • 10" 3 


0.9978 


0.89 


4.004 


2.97 • 10" 3 


0.9954 


0.28 






Shapiro test 


®Jl;Ji 


Shapiro test 


a 


mean 


sd 


W 


p-val 


mean 


sd 


W 


p-val 


2 


1.999 


5.51 • 10" 3 


0.9907 


0.72 


2.002 


6.69 • 10" 3 


0.9876 


0.25 


3 


3.001 


5.66 • 10" 3 


0.9989 


0.95 


3.001 


6.98 • 10" 3 


0.9979 


0.53 


4 


4.001 


5.59 • 10- 3 


0.9960 


0.42 


4.001 


6.40 • 10~ 3 


0.9974 


0.78 



Table 4: Sample means, variances and Shapiro Wilk Gaussianity test of cii, 
and aj L , for different values of ao 
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Full Band Estimates 



Narrow Band Estimates 



T 
2 

Fourier 



"i i 



i 



T" 



IK 

i 



-2 2 
Fourier 




r 



"T 



T 



"T 



"T 



-2-1012 
Needlet 




T 



T 



"T 



T" 



-2-1012 
Needlet 



Figure 2: Comparison among normalized distribution of Full Band and Narrow 
band estimates - a — 2, L = 1024, L\ = 724. 
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A Auxiliary Results 



This Section presents some results, mainly focussed on the behaviour of the 
Kj (a), which will be useful to develop consistency and asymptotic behaviour 
of the estimator ( 10 ). We remark that all these results hold under Condition [5J 

Proposition 24 Let Kj (a), Kji (a) and Kj 2 (a) be as in Definition^ Then 
we have: 

Kj(a) = Br°* (J s I (B,a) + o j (l) S ) ; (25) 

KjiCt0) = B^-^{I{B,a 0l a) + O] {l)), (26) 



where 



I (B,a) = 2 b 2 {u)u l - a du , (27) 
Jb- 1 

I (B, a , a) = —— — - (28) 
I (B, a) 



so that 
and 



< ci < I (B) < c 2 < +oo 



[I(a Q ,a)-i\=C I \a-a \ (29) 
Moreover, for Kj^i (a) and Kj y 2 (a) as above, we have: 

K jA (a) + Kj (a) log^ - - — B~ aj {h (B) + Oj (1)} ; (30) 

cb 

K jl2 (a) + Kj (a) log 2 B J + 2K jtl (a) log B j = —B~ aj {I 2 (B) + Oj (1)} (31) 

cb 

where 

r-B r-B 

h {B) = 2 / b 2 (u) M 1 "" (log u) du , I 2 (B) = 2 / b 2 (u) u 1 ^ (log u) 2 du, 
Jb- 1 Jb- 1 

and 

< c 3 < I 2 (B) < c 4 < +oo . 
Proof. Recalling that Nj = cbB 2j from simple calculations lead to 

1 o-a," ^ ,2 ( l\ 1 (2i + 1) l~ a 



B 3 + l 



B 



Bo / Bi Bo B- a 3 



= B-^-^-hJ b 2 (u) u^du + 0j (1) 



= B-^—{I (B)+ Oj (l)}: 
cb 
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by applying the Lagrange Mean Value Theorem we obtain 

I a {B,a) = 2b 2 (i)e- a (B-B- 1 ) forCe [B' 1 , B] 
which is a non-zero, finite positive real number. Obviously: 



Kj (a) 



^■{/o(B,a) + 0j -(l)} = Bj(ao _ a) 

TD — f\ r\ ') If /D n . \ I « /111 L V ? J 



Jf,(ao) B-^U (B,a )+ Oj (l)} 
Because by construction < C min < b 2 (£) < C max < oo, we have 
C^uB 1 - 1 * (B - B- 1 ) < I (B,a) < C^B 1 '" (B - B' 1 ) 

so that 



a )+ Oj (l)} 



max j^ a —g 



<I(B,a ,a)< ^-B 



max RQ -ct 
min 



Hence, fixing C m i n /C max < Cj < C max /C m i n , we obtain (29). We recall more- 
over that J (£?, a, a) = 1. 

As far as Kj t \ (a) is concerned, we prove (30 1. In fact: 

K jA (a) + Kj (a) logB> = 



B j+i 



cb 
B-<xj 

CB 



E * 2 



( \ 1 (2Z + 1) r a 
~E0 j RJ & B~ a i 



log 



Bi 



B 



B- 1 



2b 2 (u) u 1 - (log u) du + oj (1) 



1 



= -B-» { h (B)+o j (I)} . 
cb 

Now, we have, by applying again Lagrange Mean Value Theorem: 
h (B) = 2b 2 (0 log (0 {B- 1 - B) for £ e [B~\ B] , 
where — oo < —c% < I± (B) < C2 < +oo, so that 

1 



K Jtl (a) 



-B~ a3 — (I (B) log B> - h (B) + o, (1)) 



(32) 



Similarly, we obtain for (31 1: 



K h2 (a) + K 3 (a) log 2 B 3 + 2K j , l (a) logB j 



b^ 1 

B -^— y b 2 

C B ^ 



B' 



1 



,-B 



i \ i (2/ + i) r a 

~B~i j RJ B~3 B- a i 



log' 



BJ 



cb Jb- 1 



2b (u) u 1 -* (log uf du + oj (1) 



B-«i {I 2 ( B ) + 0J (1)} 
cb 
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which is trivially strictly positive and bounded. Hence we have: 

K j>2 (a) = B~ a i (l (B) (logB J ) 2 + 2h{B) logS^ + h (B) + 0j (1)) . (33) 

■ 

We now provide some further auxiliary results on the function Kj (a) ; these 
results are exploited in the proofs for consistency and elsewhere. 

Corollary 25 As j — > oo, we have: 

%§ 1 = - 1<,gBi + l§ = logB ' < - 1+<,,(1,,; (34) 



2 



#j\2 (a) 



^ (a) J K 3 (a) 



= (log^) 2 (l + 0j (l)) 
Proof. From we obtain: 



^-^(i + ^+.W). 

Aj (a) 7 (-BJ 

m view of 

Thus 

\ \ K 3 (a) J Ki (a) j 

■ 

Remark 26 It is immediate to see that under Condition [1| equation Q) be- 



comes ( using Proposition 24 ) 



^Nj G Q K S (a ) 
while under Condition^ we have 



) =1 + 0(B- J ) , (36) 
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Proposition 27 Let s > 0, B > 1, J x < J L Then: 



Jl 

j=.h 

Jl 

J2 B sj \ogB J 

j=Ji 



X>^' 2 log 2 B 

j=Ji 



B 



B s 



— - (b sJl -b 3 -' 1 - 1 ) 



B 



-i 1 ^ b ( bsJl { Jl -b^i) 

ss(J1 - 1) (^- 1 )-^l)) . 

L ' log 2 b(b^((j l 1 



(38) 
(39) 



B s - 1 



B h 



B s -l) (B s -1) 2 



\((J 1 - 1) 



B s 



b s -\) (B s -iy 



(40) 



Proof. The first result is trivial: 



£ s - 1 B s -I 



£ //- = £ //-< - £ /;-•• ! ! ' ~ 

j=Ji 3=0 j=0 

B J_ ^ B SJ L _ BS (J!-1)^ 



Likewise, we obtain: 



d B s _-B s \ogB 
dsB s - 1 ~~ (£« _ x) 2 



^ B sj logB j 
j=Ji 



d 
ds 



Jl 



^2 ex P{ s J lo S B } 



J=Jl 



d 

ds 



B s 



B s 



-A 



B s Jl _ ^s(Ji-l) 



logB ( B sJ - (j L —) - £? s ( Jl_1 ) ( (J, - 1) — 

B s - 1 6 V V S s - I J V 7 B s - 1 



Finally, we have 

Jl 



£ B'if ilogBf = {^{B^ - B^}\ 

j=Ji 
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-B s log 2 B 
(B s - l) 2 

B s log 2 B 

H 5 



sJl 



^-bH-^^-b^I 



b sJl [JI-t^Jl b 



b s - 1 (b s - iy 



B s log 2 B 



B s 



(b s - iy 



2 



(5 s - I) 2 



The next result combines (|38f, (|39j), (|40j) 
Corollary 28 Lef 



Jl 



Jl 



Zj^ («) = E ^ E 5 SJ j 2 log 2 B - E ^'ilogS 



XTien we have: 



Z Jl;Ji ( s ) = 



B s log B V 
B s - 1 J 



(B« - l) 2 



B 8 ^ - B^' 71 - 1 )) - B'^+Ji- 1 ) ( J L - (J x - l)) 2 

(41) 



Moreover if Ji = 1 , we ftave 

z, 7i ( S )= ( E^j (e^'/"'''^'') - (E 5 ^' lo g 5 



that 



B 3s 



lim B~ zsjL Zj r (s) = \og z B,„ 



Ji— s-oo 



(42) 



Proof. Recalling (38 1, (39) and (40), we have 



X>M (e /'-'.r i'^/' 
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B s log B 



B s - 



"\ (B' J *- - B'V*- 1 )} x Ib sJ - (j L 



B s 



B S -1J (B s -iy 



B s(.h-1) ( ( J : - 1) - 



1 



B s 



B s log B 

B s - 1 



B s -lJ (B'-iy 

2 



B 2sJ L j 



1 



B s 



+S 2s(Jl - 1) (((Ji-l)-^ L I 



b s -ij (B s -iy 

2 B s \ 



{b s - iy 



while, on the other hand: 

J2B sj j log B 



2B S 



(b s - iy 



J =1 



B s log B 

B s - 1 



s s - 1 



B s - 1 



2S «(^+Ji-i) ( J L - -J—) ( (Jl _ l) - ' 



B s -1 



B a - 1 



so that: 



Clearly if J x = 1 



B 



(b sJl - b sJi ) - b s ( Jl+j ^ (j l - jyy 



B 3s 



Z Jl (j,) = B 2s ^ log 2 g (gfl _ 1)4 + o (S 2sJi ) , 



as claimed. 



Lemma 29 Lei ylj and Bj be defined as in (19) and (20). As Jl — > oo 



\ ^1 I2 i>3 1*4 ) 



O, 



Jl \ B 2Jl 
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Proof. It is readily checked (see also [13]) that 

cum\c u C u C u Ci}=O(l- 3 r ia0 ) . 

Let us compute: 

h ' h J3J4 ' 1 N h G Q K h (a ) ' N h G K h (a ) ' N h G K h (a ) ' N j4 G K j4 (n () ) 



^ Nji G K ji{ a )) CUm (? E ? E * 



2 



- e ln TO )°"^ AiU - 

hMM.h \»=1 J ' J * v y / 

= fe |[[B (l, »- 2te i 2 (i)) B(2 ~ 4a ° )3 (l 1 ^")^ 
= o( B -vf[5f) . 



Then we have 



\GjJjxq) Gj L (a ) Gj L (a ) Gj L (a ) \ 
I Gq Go Go Go J 



\ 31323331 J 

As in |13| . the proof can be divided into 5 cases, corresponding respectively to 

l_ 31 32 33 34 ) \ jl 32 33 34, 

-^UZ cum \ J2 A n > J2 B n . Z B *> , Z B ^ [ ' ^IJZ cum { Z ^ ' J2 A h ' Z B ^ > Z B . 

^ jl J2 J3 .74 J [ jl 32 33 34 
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and 

1 



B 4Jl 



cum 



1 32 33 34 



where we have used 19 20 We have for instance 



O 



Ij^cum | ^2 A h , ^ A h > E A k ' E A ^ \ 

l_ 31 32 33 34 ) 

( u l E II ^"')<1,,,) 

\ 31,3233,34 »=1 / 

(sk E lo S 4 5 ^- 6j j =°(^E l0 S 4 = 0( 



= o 

and 

E />> .' ' E />> " ' E 11 ■ ' E B 



log 4 B Ji 
B 2Jl 



1 

cum 



J4 

Jl J2 J3 J4 

4 ^ 

Gj!,(ao) Gjjao) Gj l (q ) Gj^ao) 



CUTO ■ 



-B 4Jl I . I Go Go Go Go 

^ 3 

= O (log 4 B Jl B~ 2Jl ) ; 
The proof for the remaining terms is entirely analogous, and hence omitted. 
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